Conductivity of graphene with resonant and non-resonant adsorbates 
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We propose a unified description of transport in graphene with resonant and non-resonant ad- 
sorbates. This allows to test the validity of many approximate theories and in particular we show 
how semi-classical approximations fail near the Dirac energy for resonant and non-resonant adsor- 
bates. Close to the Dirac energy an empirical formula for the conductivity, valid in a large range 
of values of the inelastic mean- free path Li, is obtained. This formula fully takes into account lo- 
calization effects and consequences for the temperature dependence of the conductivity and for the 
magneto-conductivity are detailed. 
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Electronic transport in graphene [3-|H is sensitive 
to static defects that are for example frozen ripples, 
screened charged impurities, or local defects like va- 
cancies or adsorbates [5|-|8(. Adsorbates, which can be 
organic groups or adatoms attached to the surface of 
graphene, are of particular interest in the context of 
functionalization which aims at controlling the electronic 
properties by attaching atoms or molecules to graphene 
|9h12||. Therefore developing a correct theory of conduc- 
tivity in the presence of such defects is of major inter- 
est. Theoretical studies of transport in the presence of 
local defects have dealt mainly with semi-classical ap- 
proaches either with the Boltzmann formalism or with 
self-consistent approximations 0, [l3l - l20j . These ap- 
proaches indeed explain some experimental observations 
such as the quasilinear variation of conductivity with con- 
centration of charge carriers 9 12] or a non zero conduc- 
tivity of the order of e 2 /h at the Dirac point. 

Yet these theories have important limitations and can- 
not for example describe the localization phenomena that 
has been reported in some experiments 0, 0] • Indeed in 
the presence of a short range potential, such as that pro- 
duced by local defects the electronic states are localized 
on a length scale £ 2jJ, |22j . Therefore a sample will be 



insulating unless some source of scattering, like electron- 
electron or electron-phonon interaction, leads to a loss of 
the phase coherence on a length scale Li < £. The inelas- 
tic mean-free path Li therefore plays a fundamental role 
in the theory of conductivity of graphene with adorbates. 

In this letter we develop a numerical approach for the 
conductivity of resonant and non-resonant adsorbates on 
graphene. This method takes fully into account the effect 
of Anderson localization and allows to compute the con- 
ductivity as a function of the inelastic mean-free path 
Li. Therefore the validity of standard approaches like 
the semi-classical Bloch-Boltzmann theory and the self- 
consistent theories can be tested. In particu- 
lar for resonant scatterers we conclude that the available 
mean-field theories fail to reproduce the physics in the 



region of increased density of states close to the Dirac 
point. When the inelastic mean free path Li and the 
elastic mean-free path L e are identical we find that the 
conductivity takes a universal value a ~ 2Gq/tt in some 
energy range, close to the Dirac point where Go = 2e 2 /h 
is the quantum of conductance. For a large range of 
values of Li with Li > L e , the conductivity is given 
by a ~ (2G /it) — (aGo)Log(Lj/ L e ) where a depends 
on the type of adsorbate. From this result one can 
discuss the temperature dependence of the conductivity 
and the magneto-conductivity since the magnetic length 
L(B) — y/h/eB plays the role of a finite coherence length 
Li. 

Model for resonant and non-resonant adsorbates 

The carbon atoms of the graphene sheet that bind with 
the adsorbates create a covalent bonding and weakly con- 
tribute to the 7r-band of itinerant electrons. Therefore a 
generic model of adsorbates is obtained by removing the 
p z o rbitals of some carbon atoms of the graphene sheet. 
llMIHSl- For example an hydrogen adsorbate can 



be modeled by removing the p z orbital of the carbon atom 
that is just below the hydrogen atom. This is the model 
of resonant adsorbate that we consider here. The non- 
resonant model is constituted by two neighboring missing 
orbitals (di- vacancy). Indeed in that case the balance be- 
tween the two sub lattices of the honeycomb lattice is pre- 
served and therefore no resonant state is produced. This 
model represents a larger adsorbate that binds to the two 
neighboring carbon atoms. Let us emphasize that these 
models of missing p z orbitals do not represent real va- 
cancies. Indeed real vacancies induce large local atomic 



relaxations [12|, |25| which are not modeled here. Finally 
we consider here that the up and down spin are degen- 
erate i.e. we deal with a paramagnetic state. Indeed the 
existence of a magnetic state for various adsorbates, like 
hydrogen for example, is still debated j26[. Let us em- 
phasize that in the case of a magnetic state the up and 
down spin give two different contributions to the conduc- 
tivity but the individual contribution of each spin can be 
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analyzed from the results discussed here. With these as- 
sumptions the generic model Hamiltonian for adsorbates 
writes: 



H 



(id) 



(1) 



where represents nearest neighbours pairs of occu- 
pied sites and t — 2.7 eV determines the energy scale. In 
our calculations the vacant sites (resonant adsorbates) 
or the di-vacant sites (non-resonant adsorbates) are dis- 
tributed at random with a finite concentration. 
Evaluation of the Kubo- Greenwood conductivity 
The present study relies upon the evaluation of the 
Kubo-Greenwood conductivity using the Einstein rela- 
tion between the conductivity and the quantum diffusion. 
We evaluate numerically the quantum diffusion X 2 (E,t) 
of states of energy E for the Hamiltonian (JTJ. X 2 (E,t) 
is defined by: 



X 2 (E,t) = ((X(t) 



X{Qf) E 



(2) 



where X(t) is the position operator along the re- 
direction in the Heisenberg representation and ( ) e means 
an average over states with energy E, usually the Fermi 
energy. X 2 (E,t) is numerically evaluated using the 
MKRT approach [27r|3ljj. This method allows very ef- 
ficient numerical calculations by recursion in real-space. 
It has been used to study quantum transport in disor- 
dered graphene, chemically doped graphene, graphene 
with functionalization and graphene with structural de- 
fects 111, [l2|, |25j, I32M3 7I | . Our calculations are performed 
on samples containing up to 10 8 atoms which corresponds 
to a typical size of about one micron square. This allows 
to study systems with characteristic elastic and inelas- 
tic lengths of the order of a few hundreds nanometers. 
With characteristic lengths of such size it is possible to 
treat systems with low concentrations of adsorbates that 
are of 0.1%, 0.2%, 0.4% for resonant adsorbates (mono- 
vacancies) and of 0.5%, 1%, 2% for non-resonant adsor- 
bates (di- vacancies) . For the results presented here the 
energy resolution is of the order of 10~ 2 eV. 

The effect of the inelastic scattering is treated in a 
phenomenological way. We introduce an inelastic scat- 
tering time Tj, beyond which the propagation becomes 
diffusive due to the destruction of coherence by inelastic 
processes. This relaxation time approximation (RTA) 
has been used succesfully to compute [H| conductivity 
in approximants of quasicrystal where quantum diffusion 
and localisation effect play a essential role (39l - l4l| . Fol- 
lowing previous works [H, El] , we assume that the 
velocity correlation function Ci (E, t) of the system with 
inelastic scattering is given by 



d(E,t) ~ C(E,t)e 



-|*IA 



(3) 



where C(E, t) is the velocity correlation of the system 
without inelastic scattering. As shown in Refs. HI, 42, 43 1 
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FIG. 1: Densities of states as a function of energy, (a) 
For resonant adsorbates (mono-vacancies) with concentra- 
tions (empty circles) 0.1%, (empty square) 0.2%, (filled cir- 
cle) 0.4%. (b) For non-resonant adsorbates (di-vacancies) 
with concentrations (empty circles) 0.5%, (empty squares) 
1%, (filled circle) 2%. (Dashed lines) graphene without ad- 
sorbate. The inset shows the electron density per atom n e as 
a function f the energy with the same symbols and colors. 



the propagation given by this formalism is unaffected by 
inelastic scattering at short times (t < n) and diffusive 
at long times ( t > r») as it must be. We consider also 
that the inelastic scattering does not affect the density 
of states n{E F ). Therefore we get: 




a(E F ,n) 
D(E F , Ti ) 

Li{E F , n ) 2 



e 2 n(E F )D{E F ,n) 
Lj{E F ,n) 2 

n 

J™X 2 (E F ,t)e- t / r >dt 



(4) 
(5) 

(6) 



where E F is the Fermi energy, n(E F ) the density of 
states, Tj the inelastic scattering time and Li(E Fl Ti) the 
inelastic mean-free path. Let us emphasize that in the 
equation ([6]) X 2 (E,t) is calculated for the system with 
Hamiltonian (JlJ which has only elastic scattering. 
Densities of states 

FigureQ]shows the total densities of states as a function 
of energy for the two types of adsorbates. The case of the 
mono-vacancy has already been studied in the literature 
and our results are consistent with previous ones. In 
particular the density presents a peak around the Dirac 
point which is reminiscent of the mid-gap state produced 
by just one missing orbital (resonant scatterer). For the 
di- vacancy the density of states presents no peak because 
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FIG. 2: Conductivity a(E,Ti) and inelastic mean-free path 
Li(E, n) versus inelastic scattering time n. (a) Concentration 
0.2% of resonant adsorbates (mono-vacancies) for energies 
E = (dashed line), E = 0.04 eV (thin line) and E = 0.8 eV 
(thick line), (b) Concentration 1% of non-resonant adsorbates 
(di- vacancies) for energies E = (dashed line), E — 0.1 eV 
(thin line) and E = 1.5 eV (thick line). 



there is no resonant mid-gap state for such model (non- 
resonant scatterer). 

Conductivity versus inelastic scattering time r% and mi- 
croscopic conductivity 

Figure @ presents the conductivity a(E,Ti) and the 
inelastic mean- free path Li(E, Tj) calculated as a func- 
tion of inelastic scattering time Ti for different energies 
E. Conductivities are calculated for large values of up 
to a few 10 _11 s which corresponds to inelastic mean-free 
paths up to several hundred nanometers. Let us empha- 
size again that this is possible because we can treat large 
systems containing up to 10 s atoms. We have checked 
the convergence of our calculations with respect to the 
size of the sample. 

We define the microscopic conductivity o~m as the max- 
imum value of the conductivity over all values of r, . Ac- 
cording to the renormalization theory this value is ob- 
tained when the inelastic mean- free path Li(ji) and the 
elastic mean free path L e are identical. For larger Tj one 
has Li(ri) > L e and a(E,n) decreases with increasing 
n due to quantum interferences effects. When the lo- 
calization length £ is much larger than the elastic mean- 
free path £,/L e 1 the conductivity a(E,Ti) presents a 
plateau as a function of Tj, as long as £ ^> ii(r^) > L e . 
This is clearly observed for the highest conductivities in 
Figure ((2]). According to the standard scaling theory of 
localization this is the case when cm 3> Go = 2e 2 /h. In 
this regime o~m coincides with the semi-classical estimate 
of the conductivity as we find below. 
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FIG. 3: Microscopic conductivity o m versus electron density 
per atom n e for (a) resonant adsorbates (mono-vacancies) and 
(b) non-resonant adsorbates (di- vacancies), for 3 concentra- 
tions (see caption of figure [I}. [Inset provide a zoom on the 
low concentration limit.] The predictions of the Boltzmann 
theory are represented by dashed lines. For (a) the Boltzmann 
conductivity is obtained with the self-energy of fiol j . For (b) 
we just provide the Boltzmann value close to the Dirac energy. 



Validity of semi- classical and self-consistent theories 
Figure © represents the values of the microscopic con- 
ductivity cm as a function of the number of electrons per 
carbon atom n e . The comparison with the conductivity 
ob calculated with the semi-classical Bloch-Boltzmann 
approach confirms that for large values <jm 3> Go then 

For the case of resonant adsorbates (mono- vacancies) 
around the Dirac point there is an energy region in which 
ctm is of the order of 4e 2 /wh = 2Go/n. The maximum 
energy of this region corresponds to about one electron 
per impurity. This is consistent with predictions of self- 
consistent theories [j| fl3l - [20l ] . Yet inside this plateau a 
peak of conductivity appears very close to the Dirac point 
which is not predicted by self-consistent theories. This 
peak is not obtained by [lj| and is present in the cal- 
culation of [2(J although much less marked than in the 
present work. We believe that the high energy resolu- 
tion of the present calculation (~ 10~ 2 cV) explains the 
disagreement with the two previous studies. 

For the case of non-resonant adsorbates (di-vacancies) 
the conductivity minimum appears in a very narrow en- 
ergy range. Again the proper description of this narrow 
structure requires a high energy resolution and indeed we 
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FIG. 4: Conductivity a in units of Go as a function of inelas- 
tic scattering length Li for the 3 concentrations of adsorbates 
(with the same symbols as in figure [TJ (a) resonant adsor- 
bates (mono- vacancies) at energies E = 0.03 eV (dashed line), 
and E — 0.04 eV (continuous line), (b) non-resonant adsor- 
bates (di- vacancies) at E = OeV. The dot-dashed straight 
lines show the slope —a for Li ^> L e (see text): (a) a = 0.25 
and (b) a = 0.75. The inset shows a(Li) at E = in a) and 
a(Li) at E= 0.1 eV in b). 



checked that with a lower energy resolution this plateau 
disappears 34]. This narrow minimum is absent from 



the calculations on the epoxy group in [111 ]. Again we 
attribute this difference with the present work to a lower 
energy resolution in 11 1 . 

Conductivity versus inelastic mean free path and 
magneto- conductivity 

We consider now the case when the inelastic mean- 
free path Li is greater than the elastic mean-free path 
L e (figure [4|. Let us discuss first the energies where the 
microscopic conductivity cfm is of the order of 4e 2 /irh. 
At these energies the conductivity is well represented by 



4e2 



2e\ (Li 

a x Log U 



(7) 



The coefficient is a ~ 0.25 for mono- vacancies and 
a ~ 0.75 for di-vacancies. This law is reminiscent of 
the result of the perturbation theory of 2-dimension An- 
derson localization for which a ~ 1/7T [44] . An important 
difference here is that <jm is of order Go which means that 
the regime is not perturbative close to the Dirac point. 



Our study gives a description of the temperature de- 
pendence of the conductivity provided that the temper- 
ature dependence of the inelastic mean free path L^T) 
is known. This will depend on the dominant scattering 
mechanism (usually electron-electron or electron-phonon 
scattering) and is widely discussed in the litterature [3] . 
Our study provides also a description of the magneto- 
conductivity since the magnetic length L(B) = ^Jh/eB 
plays the role of a finite coherence length such as the in- 
elastic mean- free path Li. The magneto-conductivity is 
Aa(B) ~ a(Go/2)ALog(B) and should be observed as 
soon as the magnetic length L(B) — y/h/eB is smaller 
than the inelastic mean- free path Li(T) which depends 
on the temperature T in a given sample. 

For the energies close to the Dirac point where cm 
tends to be larger than 4e 2 /Trh the conductivity does not 
follow the above law with respect to Li. In the mono- 
vacancy case a divergence of the localization length at the 
Dirac energy is predicted in 2l|. Here on the contrary 
the conductivity tends to zero at large inelastic length 
which indicates a finite localization length. This discrep- 
ancy can be explained by the finite energy resolution of 
our calculation which does not allow to capture the di- 
vergence at E = 0. Since this resolution is of the order 
of 10~ 2 cV, wc conclude that the divergence of the local- 
ization length exists only in a narrow energy range and 
could be difficult to observe experimentally. 

Conclusion 

In summary we have analyzed the conductivity of 
two generic models representative of resonant and non- 
resonant adsorbates on graphene. For two dimensional 
models with Anderson localization, like the ones studied 
here, the inelastic mean-free path Li is a key parameter 
and all the calculations here are done as a function of 
this characteristic length scale. These calculations allow 
to test the validity of many approximate theories devel- 
oped so far. In particular, for resonant absorbates, the 
available mean-field theories do not describe correctly the 
electronic structure and the physics of transport in the 
central peak of the density of states, at the Dirac energy. 

Close to the Dirac energy the microscopic conductiv- 
ity <7M) i-e. the conductivity obtained when the inelastic 
mean free path Li and the elastic mean-free path L e are 
identical, is of the order of 4e 2 /Trh. As we show this 
occurs in a region of energies which is larger for reso- 
nant than non-resonant adsorbates. When the micro- 
scopic conductivity ctm is of the order of 4e 2 /irh a linear 
variation of the conductivity with Log(Li/ L e ) is calcu- 
lated. This allows some predictions for the temperature 
dependence of the conductivity and for the magneto- 
conductivity. The prediction of a divergence of the lo- 
calization length at the Dirac energy [2l[ cannot be con- 
firmed by our calculations. This suggests that this diver- 
gence occurs on an energy range less than about 10~ 2 eV 
at the concentrations considered here. 
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